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ABSTRACT 

Using a statistically rigorous analysis method, we place limits on the ex- 
istence of an isotropic stochastic gravitational wave background using pulsar 
timing observations. We consider backgrounds whose characteristic strain spec- 
tra may be described as a power-law dependence with frequency. Such back- 
grounds include an astrophysical background produced by coalescing supermas- 
sive black-hole binary systems and cosmological backgrounds due to relic grav- 
itational waves and cosmic strings. Using the best available data, we obtain 
an upper limit on the energy density per unit logarithmic frequency interval of 
^SMBHj^j^^g y]^~)/^2 < ]^ g X 10~^ for an astrophysical background which is five times 
more stringent than the earlier Kaspi et al. (1994) limit of 1.1 x 10~^. We also 
provide limits on a background due to relic gravitational waves and cosmic strings 
of Qf'^l/S yT)h^ < 2.0 x 10"^ and nf{l/8 yi)h^ < 1.9 x 10"^ respectively. All 
of the quoted upper limits s correspond to a 0.1% false alarm rate together with 
a 95% detection rate. We discuss the physical implications of these results and 
highlight the future possibilities of the Parkes Pulsar Timing Array project. We 
find that our current results can 1) constrain the merger rate of supermassive bi- 
nary black hole systems at high red shift, 2) rule out some relationships between 
the black hole mass and the galactic halo mass, 3) constrain the rate of expan- 
sion in the inflationary era and 4) provide an upper bound on the dimensionless 
tension of a cosmic string background. 
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1. Introduction 

Pulsar timing observations (see Lorimer & Kramer 2005, Edwards, Hobbs & Manchester 
2006 for a review of the techniques) provide a unique opportunity to study low-frequency 
(10-9 - 10-^ Hz) gravitational waves (GWs) (e.g., Sazhin 1978, Detweiler 1979, Bertotti, 
Carr, & Rces 1983, Foster & Backer 1990, Kaspi, Taylor & Ryba 1994, Jenet et al. 2005). 
Sources in this low-frequency band include binary supermassive black holes, cosmic super- 
strings, and relic gravitational waves from the Big Bang (Jaffe & Backer 2003, Maggiore 
2000). 

An isotropic stochastic background can be described by its characteristic strain spectrum 
hcif) which, for most models of interest, may be written as a power-law dependence on 
frequency, /: 

/ 



Table 1 shows the expected values of A and a for different types of stochastic back- 
grounds that have been addressed in the literature. The characteristic strain is related to 
the one-sided power spectrum of the induced timing residuals, P{f), as 

Pif) = ^^Jshcif?, (2) 



and to D,gu,{f), the energy density of the background per unit logarithmic frequency interval, 
as 

^gM) = l^/h^(f)" (3) 

where Hq is the Hubble constant. Note that the one-sided power spectrum, P(/), is defined 
so that 

poo 

PU)df = (4) 



Jo 



where is the variance of the arrival time fluctuations, or timing residuals, generated by 
the presence of the GW background. Since cr^ has the physical units of s^, P(f) has the units 
of s^/Hz or s^. 

Jenet et al. (2005) developed a technique to make a definitive detection of a stochastic 
background of GWs by looking for correlations between pulsar observations. It was shown 
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that approximately 20 pulsars would need to be observed with a timing precision of ~100ns 
over a period of 5 years in order to make such a detection if the GW background is at the 
currently predicted level (Jaffe & Backer 2003, Wyithe & Loeb 2003, Enoki et al. 2004, 
Sesana et al. 2004). The Parkes Pulsar Timing Array (PPTA) project (Hobbs 2005) is 
trying to achieve these ambitious goals, but the currently available data-sets do not provide 
the required sensitivity for a detection. In this paper, we introduce a method to place an 
upper bound on the power of a specified stochastic GW background using observations of 
multiple pulsars. Full technical details of our implementation will be published in Hobbs et 
al. (2006). Here, this method is applied to data (see Section 2) from seven pulsars observed 
for the PPTA project combined with an earlier publicly available data-set. 

Upper limits have already been placed on the amplitude of any such background of 
GWs. Using eight years of observations for PSR B1855-I-09, Kaspi et al. (1994) obtained a 
hmit of Qgh'^ < 1.1 x 10"'^, where Ho = lOO/i km s^^ Mpc^^, at the 95% confidence level^ 
for the case when a = —1 (i.e. flgu, is independent of frequency). This work was continued 
by Lommcn (2002) who used 17yrs of observations to obtain Qgh"^ < 2 x 10"^. However, 
the statistical method used for both of these analyses has been criticized in the literature 
(see, for instance, Thorsett & Dewey 1996, McHugh et al. 1996, Damour & Vilenkin 2005). 
In this paper, we develop a frequentist technique, similar to that used by the LIGO science 
collaboration (Abbott et al. 2006), to place an upper bound on A, given a. The technique 
makes use of a statistic, T, defined below, which is sensitive to red noise in the pulsar timing 
residual data. Upper bounds on A are determined using T together with a specified false 
alarm rate, Vf, and detection rate, Vd- Monte-Garlo simulations are used to determine these 
probabilities by generating pulsar pulse times-of-arrival consistent with a GW background. 
All of the upper limits quoted in this paper correspond to Vf — 0.1% and Vd — 95%. 



^Thcir more stringent constraint of flgh^ < 6 x 10^^ was obtained when data from PSRs B1855+09 and 
B1937+21 were combined. Since the data from PSR B1937+21 is far from white, we believe this Hmit is 
artificially low and therefore restrict our discussion to the PSR B1855+09 data only. 



Table 1: The expected parameters for predicted stochastic backgrounds 



Model 


A 




a 


References 


Supermassive black holes 


10-15 


- 10-1^ 


-2/3 


Jaffe & Backer (2003) 










Wyithe & Loeb (2003) 










Enoki et al. (2004) 


Relic GWs 


10-1^ 


- 10-15 


-1 - -0.8 


Grishchuk (2005) 


Gosmic String 


10-16 


- 10-1^ 


-7/6 


Maggiore (2000) 
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2. Observations 

We expect the isotropic background to generate timing residuals with a "red" spectrum: 
a spectrum with excess power at low frequencies or, equivalently, long time scale correlations 
in the residuals. Therefore, we have restricted our analysis to those pulsars having formally 
white spectra: a spectrum with statistically equal power at all frequencies or no correlations 
in the residuals. This allows us to put the best upper limit on the background by bounding 
the level of any red process in those data sets. Three separate tests were used in order to 
determine the statistical properties of the data and to select data-sets that are statistically 
white. First, the normalized Lomb-Scargle periodogram was calculated for each residual 
time series. No significant peaks were seen in any of the data used. Second, the variance of 
the residuals was shown to decrease as 1/n where n is the number of adjacent time samples 
averaged together. If the data were correlated the variance would not scale as 1/n. Third, 
no significant structures were seen in the polynomial spectrum (defined below) for each 
individual spectrum or in the averaged spectrum. Note that the publicly available data-set 
for PSR B1937-I-21 (Kaspi et al. 1994) was not used in our analysis since its timing residuals 
do not pass these three tests. 

We made use of the following data-sets which passed the tests: 1) observations of 
PSR B1855+09 (also known as PSR J1857+0943) from the Arecibo radio telescope that 
are pubhcly available (Kaspi et al. 1994), 2) observations for PSRs B1855+09, J0437-4715, 
J1024-0719, J1713+0747, J1744-1134, J1909-3744 and B1937+21 (J1939+2134) using the 
Parkes radio telescope and reported by Hotan et al. (2006) and 3) recent observations of 
all of these pulsars made as part of the PPTA and related Swinburne timing projects. The 
Kaspi et al. (1994) data-set was obtained at ~1400MHz over a period of Syrs. The PPTA 
observations, which commenced in February 2004, include ~20 millisecond pulsars and use 
the 10/50-cm dual- frequency receiver and a 20-cm receiver at the Parkes radio telescope. 
Each pulsar is typically observed at all three frequencies in sessions at intervals of 2 - 3 
weeks. The results used here were obtained using a correlator with 2-bit sampling capable 
of bandwidths up to 1 GHz and a digital filterbank system with 8-bit samphng of a 256 MHz 
bandwidth. The PPTA observations and the earlier Hotan et al. (2006) data-sets also used 
the Caltech Parkes Swinburne Recorder 2 (CPSR2 - see Hotan et al. 2006), a baseband 
recorder that coherently dedisperses two observing bands of 64 MHz bandwidth, centered 
on 1341 and 1405 MHz for observations at 20 cm and around 3100 MHz and 685 MHz for 
(simultaneous) observations with the coaxial 10/50 cm receiver. Full details of the PPTA 
project will be presented in a forthcoming paper; up-to-date information can be obtained 
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from our web-site^. Unfortunately, our stringent requirements on the "whiteness" of the 
timing residuals has restricted the use of some of our nominally best-timing pulsars. For 
instance, even though a 10-yr data span is available for PSR J0437— 4715, the full-length 
data-set is significantly affected by calibration and hardware-induced artifacts as well as 
other unknown sources of timing noise. 

A listing of the pulsars observed, the observation span, number of points and weighted 
rms timing residual after fitting for the pulsars' pulse frequency and its first derivative, 
astrometric and binary parameters are presented in Table 2. Arbitrary offsets have been 
subtracted between data-sets obtained with different instrumentation. Combining these 
data-sets provides us with data spans of ~20 yr for PSR B 1855-1-09 and ~ 2 — 4 years for the 
remaining pulsars. The final timing residuals are plotted in Fig. 1. 

3. New upper bounds on the stochastic backgrounds 

The goal here is to use the measured timing residuals from multiple pulsars in order to 
determine the smallest value of A that can be detected for a given a as defined by equation 1. 
This is done in a three-step process. First, a detection algorithm is defined that is sensitive 
to the presence of the background. Second, this algorithm is tuned so that in the absence of a 
signal (i.e. A — 0), the probability of the detection scheme falsely detecting the background 
is set at Vf, known as the false alarm rate. Lastly, for the given detection scheme and 
false alarm rate, the upper bound, ^upper, is chosen so that the probability of detecting a 



^http : / /www. atnf . csiro.au /research /pulsar / ppta 
Table 2: Pulsar observations used for this analysis 



Pulsar 


Telescope 


Span 


N 


rms residual 






(d) 






J0437-4715 


Parkes 


815 


233 


0.12 


J1024-0719 


Parkes 


861 


92 


1.10 


J1713+0747 


Parkes 


1156 


168 


0.23 


J1744-1134 


Parkes 


1198 


101 


0.52 


J1857+0943 


Arecibo/Parkes 


7410 


398 


1.12 


J1909-3744 


Parkes 


866 


2859 


0.29 


J1939+2134 


Parkes 


862 


231 


0.21 
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Fig. 1. — Pulsar timing residuals. The length of the vertical hne on the left hand edge 
represents 10 /is. 
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background with A — ^upper is Vd- For this paper, the false alarm rate is set to 0.1% while 
the upper-bound detection rate is set to 95%. 

Since all current models of the background predict that the induced timing residuals 
will be "red" (the spectrum increases at lower frequencies), the detection scheme employed 
here is defined to be sensitive to a red spectrum. The existence of a red spectrum in the 
timing residuals is therefore necessary, but not sufficient, evidence for the existence of a GW 
background. Hence, we can use a statistic sensitive to a red spectrum in order to place 
an upper bound on the amplitude of the characteristic strain spectrum. Since the data- 
sets are irregularly sampled and cover different time spans, a spectrum based on orthogonal 
polynomials is employed. Each pulsar data-set consists of Up measured residuals, Xp{i), a 
time tag tp{i), and an uncertainty (Jp{i), where i is the data sample index and p is an index 
referring to a particular pulsar. The time tags are scaled so that normalized time tags, Tp{i), 
run from —1 to 1. These Tp{i) values are used in a weighted Gram-Schmidt orthogonalization 
procedure to determine a set of orthonormal polynomials defined from 

where is the Tth order polynomial evaluated at Tp{i) and is the standard Kronecker 
delta function. Note the highest power of t in is I. For the case where r is continuous 
and (7p{i) — 1, the above sum becomes an integral and jp{i) become the familiar Legendre 
polynomials. The following coefficients are calculated using the orthonormal polynomials, 
and the timing residuals Xp{i): 

The pulsar average polynomial spectrum is given by 

P, = ^(3)! (7) 

Up 

where the weighted variance, Vp, is defined as ^ Yl7=o ~ ^pY/'^pi'^) ^ mean 

of X. Since the stochastic background is red. Pi will be large for low values of I if the 
background significantly influenced the residuals. Hence, T = I^;=o Pi be used as a 
statistic to detect the background. An upper limit of seven is used since 95% of the power 
is contained in the first seven polynomials for the case of a = —2/3. The background will 
be "detected" if T > Tq where Tq is set so that the false-alarm rate is given by Vf. 



-8- 



A Monte-Carlo simulation was used to determine Tq and ^upper- Complete details of the 
simulation and its implementation may be found in Hobbs et al. (2006), but a brief overview 
is given here. The simulation, undertaken in the pulsar timing package TEMP02 (Hobbs, 
Edwards & Manchester 2006), generates an ideal TOA data-set (with the same sampling as 
the observed data) from a measured set of TOAs and a given timing model. The fluctuations 
due to the GW background for a given A and a are introduced into the TOAs by adding 
together 10,000 sinusoidal GWs which come from random directions on the sky and have 
randomly chosen frequencies in the range (l/2000yr)-(l/0.5d). As a test of the simulation, 
the ensemble- averaged power spectrum of the simulated residuals was calculated over a time 
scale much larger then the longest GW timescale (i.e. 2000 years) and was shown to be 
consistent with equation 2 as expected. The GW residuals are then added to the ideal TOA 
data-set for each pulsar. In order to include the effects of measurement noise, the measured 
timing residuals are added back into the data-set, but randomly shuffled. This ensures that 
the added noise has the same probability distribution as the actual measurement noise. In 
this way, a new set of TOAs are generated that include both measurement noise and the 
GW background. Note that the shuffling procedure is only valid when the data have a white 
spectrum. Otherwise, the spectral properties of the original data set and the shuffled data 
set will not be the same. This simulated TOA data set will then be analyzed in exactly the 
same way as a real data set. Hence, all the systemic effects which inhibit gravitational wave 
detection such as low order polynomial removal. Earth's orbital motion, annual parallax 
effects, and orbital companion effects are appropriately accounted for in the simulation. 

To calculate Tq, the simulation generates 10000 independent simulated sets of TOAs 
for each pulsar with A = (i.e. no GW background). The statistic T is calculated for each 
of the 10000 trials. Using this set of T values together with the chosen false alarm rate, Vf, 
the value of Tq can be determined. Once Tq is chosen, the simulation is used to generate 
TOA data-sets that include the effects of GWs. For a given value of A, the probability of 
detection is determined using T and Tq. ^upper is chosen to be that value of A when the 
probability of detection is equal to Vd- 

Note that the effects of unknown time offsets ("jumps") in the data-sets are included in 
the calculation of both Tq and A^pper using this technique since tempo2 fits for these offsets 
in the TOA data-set after the GW background has been added. Since we are using TEMP02 
to analyze the data, the effects of all the fitting procedures are being taken into account. 
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1-2 -1 1 

a 

Fig. 2. — Minimum detectable A (or flgyj(lyT~^)h'^ - right axis) versus a for our current 
limits (solid line) and potential future limits from the PPTA (dashed line). The star symbol 
indicates the limit obtainable using the Kaspi et al. (1994) observations of PSR B1855+09. 
From left to right the near-vertical dotted lines indicate the expected range of amplitudes 
for the cosmic strings, relic GW and supermassive black hole background respectively. 
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4. Results 

Using the pulsar data-sets described above, the 95% detection rate upper bound with a 
false alarm rate of 0.1% is given in Table 3 for different values of a. The relationship between 
A and a is shown in Fig. 2. These upper bounds on A can be converted to an upper bound 
on the normalized GW energy-density per unit logarithmic frequency interval, flgwif ) using 
equations 1 and 3. Our limits on Qg^(lyr~^) are indicated on the right-hand axis of Fig. 2. 

We can compare our results to the previously published limit of Kaspi et al. (1994) who 
obtained that ^^^(l/S yr)/i^ < 1.1 x 10"^ (star symbol in Fig. 2). Using the same data-set 
as Kaspi et al. (1994) our method provides a similar limit of ^2^^(1/8 yr)/i^ < 1.3 x 10~^. 
Combining this data-set with our more recent observations provides a more stringent limit 
of Qg^(l/8 yr)/i2 < 1.9 x 10"^ 

The most stringent hmit reported to date was obtained by Lommen (2002). Unfor- 
tunately, these observations are not publicly available. In order to compare our technique 
we use the original PSR B1855-I-09, Kaspi et al. (1994) data-set along with two simulated 
white data-sets that realistically model the NRAO 140-foot and Arecibo observations which 
form the remainder of the Lommen (2002) data (we simulate 60 observations with an rms 
residual of 5fis between MJDs 47800 and 51360 for the 140-foot telescope and a further 60 
observations with an rms residual of l^s between MJDs 50783 and 52609 for the most recent 
Arecibo data). As we simulate neither systematic effects nor timing noise our limit will be 
more stringent than could be obtainable using the real data-set. For a = —2/3, we obtained 
that A < 9 X 10~^^ corresponding to f2g^(l/17 yr)/i^ = 8 x 10~^. This limit is a factor of 
four less stringent than that reported by Lommen (2002). 

Using simulated data, the upper bounds that can be expected from future experiments 
can be determined. The goal of the PPTA is to time 20 pulsars with an rms timing residual 
of 100 ns over 5 years. The dashed hne in Fig. 2 plots A versus a for such a data-set which 
could potentially provide a limit on a background of supermassive black hole systems of 
Aupper < 6.5 X 10-1^ or Qg^il/S YY)h? < 6.6 x 10"^^ (see Table 3). 

In Jenet et al. (2005), techniques to use an array of pulsars to detect a stochastic 
background of GWs with a — —2/3 were developed.^ Given a value for A^pper, one can use 
such techniques to determine the probabihty of definitively detecting the GW background 
using the completed PPTA data-sets (20 pulsars with an rms timing residual of 100 ns over 
5 years) if A were equal to ^upper- In terms of the parameter, defined in Jenet et al. 



^Note that A as defined here is larger by a factor of V3 compared to the definition of A used in Jenet et 
al. (2005). The definition used here is consistent with Jaffe & Backer (2003) and Wyithe & Lobe (2003). 
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(2005), a significant detection would occur if 5" > 3.1. This corresponds to a 0.001 false 
alarm rate. For the case of a = —2/3, the expected value of S (assuming ideal whitening) is 
about 4.1 for A = A^pp^j.. Since the probability distribution of S is approximately Gaussian, 
the probability of > 3.1 when (S) = 4.1 is 85%. Hence, the GW background would 
be detected 85% of the time. For the case of 10 years of observations, the detection rate 
increases to over 99% of the time. 



5. Implications and Discussion 

The upper bound on the stochastic background can be used to probe several aspects of 
the Universe. Precisely what is being constrained depends on the physics of the particular 
background in question. Here, both the measured upper bounds using the currently available 
data and the expected upper bounds using the full five-year PPTA data-set are discussed in 
the context of several GW backgrounds. 



5.1. Supermassive black holes 

A GW background generated by an ensemble of supermassive black holes distributed 
throughout the universe has been investigated by several authors (Jaffe & Backer 2003, 
Wyithe & Loeb 2003, Enoki et al. 2004). In general, the characteristic strain spectrum for 
this background can be written as: 

M/,^.«o-..(^)-((^)-)"^(^)-.. ,3) 

where 

N{z) a{z) dz 
No °d(^) (1 + ^)4/3' 

a(z) is the cosmological scale factor written in terms of red shift, z, a{z) is the derivative 
of a{z) with respect to cosmic time, Hq is the Hubble constant, the chirp mass Mc = 
[MiM2(Mi + M2y^/^f^^ of a given binary system, () represents ensemble averaging over 
all the systems generating the background, N{z) is the galaxy merger rate as a function 
of red shift, and Nq is the present day number density of merged galaxies that created 
a black-hole binary system. The values of each of these physical quantities are currently 
poorly constrained and each investigator has chosen a different parameterization. Under 
the framework described by Jaffe and Backer (2003), (Mc^^) and A^o are constrained by 
observations at the current epoch to be (Mc) 2.3 x IO'^'Mq and Nq ^ 1/Mpc^. They 
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parameterized the galaxy merger rate such that R{z) goes as (1 + z^. Hence, / depends 
on 7. Combining the estimates of (Mc) and Nq together with our measured upper bound of 
^uppcr = 1.1 X 10~^^, one finds that / < 3. Using the full PPTA after 5 years, one expects 
/ < 0.8. These constraints together with the calculations of Jaffe & Backer (2003) (see 
Fig. 4 in their work) constrain 7. Currently the limit on 7 is < 2.6 and with the full PPTA 
7 < 0.4. This value is expected to lie somewhere between —0.4 and 2.3 (Carlberg et al. 2000, 
Patton et al. 2002). Current PPTA sensitivity (i.e. using the data presented in this paper) 
is just above the expected range, while the full PPTA should be able to place meaningful 
constraints on this exponent. 

In the Wyithe & Loeb (2003) work, both (Mc ) and / depend strongly on the black- 
hole versus galactic-halo mass (Mbh — A^hm) relationship. They discuss several different 
scenarios which yield different Mbh — Mum relationships and hence different levels of the 
background. For the case of an Mbh — Mhm relationship determined by Ferrarese (2002), 
the expected value of A is 2 x 10~^^. For the Mbh ~ Mhm relationship derived from Navarro, 
Frenk, & White (1997), A = 5 x 10^^^. Using an Mbh — Mhm relationship derived from 
simple considerations of BH growth by feedback from quasar activity (Wyithe & Loeb 2003, 
Haehnelt, Natarajan & Rees 1998, Silk & Rees 1998), A ^ 10^^^. Our measured upper limit 
for a = —2/3 cannot rule out any of these models. However, if only a limit is obtained from 
the full PPTA observations, it will rule out all of the models described above. 



A relic gravitational wave background is generated by the interaction between the large- 
scale dynamic cosmological metric and quantum fluctuations of the metric perturbations 
occurring in the early universe (Grishchuk 2005). In the nano-Hz frequency regime, the 
background takes the following form: 



where hc{fh) is the magnitude of the characteristic strain spectrum at / = Hq, an is the 
current value of the cosmological expansion factor, and 02 is the value of the expansion 
factor at the start of the matter-dominated epoch. Note that this expression is not valid 
in the ultra- low frequency regime where f ^ Hq. The notation used here is consistent 
with Grishchuk (2005) except for a which is related to Grishchuk's parameter j3 hj a — 
1 -\- (3. The exponent determines the evolution of the inflationary epoch which starts the 
GW ampliflcation process. When a — —1, the scale factor grows exponentially with global 
cosmic time. The ratio a^jau is believed to be about 10~^. hc{fh) is constrained by cosmic 



5.2. Relic gravitational waves 




(10) 
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microwave background measurements to be about 10~^. Using these values and assuming 
the vahdity of the amphfication scenario described in Grishchuk (2005), the upper bound 
on A may be used to constrain a. The upper bound on a is given by the solution to the 
following equation: 

hM)rA^T(^Y = A(a). (11) 



fh J V 0.// , 

The above equation yields a < —0.7 for the current PPTA and a < —0.84 for the full 
PPTA. Within the theoretical framework described by Grishchuk (2005), if a is larger than 
—0.80, small scale gravitational waves will effect primordial nucleosynthesis, while an a less 
than —1.0 will result in an infinitely large energy density in small scale gravitational waves. 
Hence, the full PPTA will be able to place useful constraints on the rehc gravitational wave 
background. Since a determines the rate of expansion in the inflationary epoch, it turns out 
that it is related to the equation of state of the "matter" in that epoch by 

- = » = (12) 

where P is the pressure and e the energy density. The full PPTA will constrain w in the 
early universe to be greater than —1.17. This would limit inflationary models based on 
"quintessence" and "phantom energy" (Nojiri et al. 2006; Padmanabhan 2005). 



5.3. Cosmic strings 

It has been proposed that oscillating cosmic string loops will produce gravitational wave 
radiation (Vilenkin 1981). Recently, Damour & Vilenkin (2005) discussed the possibihty of 
generating a stochastic GW background using a network of cosmic superstrings. Using a 
semi-analytical approach, they derived the following characteristic strain spectrum valid in 
the pulsar timing frequency range (see their equation 4.8): 

Kif) = 1.6 X 10-^V'Y''%f(h/.65y"> (^) (-^) (13) 

where /x is the string tension, G is Newton's constant, c is the average number of cusps per 
loop oscillation, p is the reconnection probability, e^ff is a loop length scale factor, and h 
is the Hubble constant in units of 100 km/s/Mpc. Note that for the above estimate, h was 
evaluated at 0.65 in order to be consistent with Damour & Vilenkin (2005). The combination 
G/i is the dimensionless string tension which characterizes the gravitational interaction of the 
strings. The predicted string tensions are 10~^^ < G/j, < 10"^ (Damour & Vilenkin 2005). 
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Using the above spectrum together with the measured upper bound on he for a — —7/6, an 
upper bound may be placed on the dimensionless string tension: 

G^i < 1.5 X W-'c-'/y/^elff{h/.65)-'/^ (14) 

As emphasized by Damour & Vilenkin (2005), the above expression for the upperbound may 
be simphfied using the fact that both p and eg// are less then one and h is expected to be 
greater than .65: 

G/i < 1.5 X 10-^3/2. (15) 

Using standard model assumption where c = 1, the upperbound becomes G/j, < 1.5 x 10^^. 
This is already limiting the parameter space of the cosmic string model of Sarangi and Tye 
(2002). With the full PPTA, the limit will become G/x < 5.36 x lO^^^ Hence, the full PPTA 
will either detect GWs from cosmic strings or rule out most current models. 
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Table 3: Current and potential future limits on the stochastic gravitational- wave background 



a 


A 




0,^(1/1 yr)/i2 


1^3^(1/8 yr)/i2 


fi,^(l/20 yr)/i2 


-2/3 


1.1 X 10" 


-14 


7.6 X IQ-s 


1.9 X 10-^ 


1.0 X 10"^ 


-1 


5.7 X 10" 


-15 


2.0 X 10-^ 


2.0 X 10-s 


2.0 X 10"^ 


-7/6 


3.9 X 10- 


-15 


9.6 X lO-'^ 


1.9 X 10-^ 


2.6 X 10-^ 


-2/3 


6.5 X 10- 


-16 


2.7 X 10-^° 


6.6 X 10-^^ 


3.6 X 10-^1 


-1 


3.8 X 10- 


-16 


9.1 X 10-^1 


9.1 X 10-^^ 


9.1 X 10-^1 


-7/6 


2.8 X 10- 


-16 


4.9 X 10-^1 


9.9 X 10-^^ 


1.3 X 10-^° 



Note. — The upper half of the table gives limits derived from current observations. Limits based on 
timing 20 pulsars with an rms timing residual of 100 ns over 5yr are given in the lower half of the table. 



